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Abstract. We study the effect of liquidity freezes on an economic agent optimizing her utility of 
consumption in a perturbed Black-Scholes-Merton model. The single risky asset follows a geometric 
Brownian motion but is subject to liquidity shocks, during which no trading is possible and stock 
dynamics are modified. The liquidity regime is governed by a two-state Markov chain. We derive 
the asymptotic effect of such freezes on optimal consumption and investment schedules in the two 
cases of (i) small probability of liquidity shock; (ii) fast-scale liquidity regime switching. Explicit 
formulas are obtained for logarithmic and hyperbolic utility maximizers on infinite horizon. We 
also derive the corresponding loss in utility and compare with a recent related finite-horizon model 
of Diesinger, Kraft and Seifried (2010). 



1. Introduction 



The theory of optimal investment has been one of the cornerstones of the mod ern mathemat ical 
finance. Working in the context of the Black-Scholes model of a single risky asset, MertonI ( 1971 ) in 
his seminal work obtained closed-form expressions for optimal investment and consumption strate- 
gies of an agent with power-utility preferences. These results relied on the fundamental assumptions 
of frictionless and continuous trading. In the past thirty five years, extensive efforts have been de- 
voted to reco nciling this assumption with real financial markets. A large body of literature (see e.g. 
the textbook iKaratzas and Shrevd (|l998l ) and references therein) treats optimal investment with 



(either proportional or fixed) tran saction costs; another significant strand addresses optimal i nvest' 



ment with discrete trading times (IRogers 



Matsumoto 



20061 : 



Cretarola et al 



201ol ). Fin ally, several authors considered optiinal investment 



2001 



Rogers and Zanel . l2002l : IPham and Tankovl . l2008l : 



2000 



Bauerle and Rieder 



20041 1 



consumption models unde r stoch astic volatility (jFleming and Hernandez-Hernandezl . l2003l : lFouque et al 



Recently, t he issue of illiquid financ i al markets has be en considered. From one direction, sev- 



eral authors (j(^etin and Rogersl . 



20071 : l(^etin et al.l . |2004| ) have analyzed optimal investment by 



a large trader whose actions affect the dynamics of the market. From a different perspective, 
Diesinger et al.l (|2010l ) have treated markets where trading may be interrupted. The latter phe- 
nomenon means that during times of crisis markets freeze and often no trading is possible at all, 
no matter the price. Thus, during major financial dislocations, trading halts and circuit-breakers 
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are implemented to alleviate panic. These shutdowns typically follow a major market drop and 
are often triggered by the index losing a preset percent of its value within a single trading ses- 
sion. For instance, on October 28, 1997 the New York Stock Exchange (NYSE) experienced two 
circuit-breakers once the Dow Jones index lost 5% and then 7.5% percent of its value. In each 
case, trading was halted for an hour. More recently, during the 2008 financial crisis, the Russian 
RSX exchange was closed by government decree for three days during Sep 20-23, 2008. Many other 
emerging market bourses in Brazil, India, Thailand, Bangladesh, etc., also faced repeated halts for 
parts or all of a trading session in September-October 2008. Similarly, after the 9/11 attacks in 
2001, the NYSE was closed for four days. On many of the smaller national exchanges, investors 
must contend with regular trading halts due to political, regulatory and financial upheavals. 

More lengthy interruptions have been witnessed during times of war. In the 1940s, the Zurich 
stock exchange closed for over a month in May 1940; the Frankfurt stock exchange was closed for 
5 months in 1945. At the onset of World War I the NYSE was closed from July 31, 1914 until 
November 28, 1914 (jPiesinger et al.l . l20ld ). Moving away from stocks, lengthy liquidity freezes are 
commonplace in more exotic asset classes. The recent credit crisis is a case in point; 2007-2008 saw 
major liquidity crises in US sub-prime mortgage backed securities (MBS), Canadian asset backed 
commercial paper (ABCP) and worldwide collaterized debt obligations (CDOs). The collapse of 
Lehman Brothers in September 2008 led to several months of virtual shutdown in trading of credit 
products and there were many anectodal stories of investors unable to sell their "toxic" assets at 
any price. 

The effect of such liquidity crises on optimal investment is threefold. First, inability to rebalance 
the portfolio during the interruption reduces the long-run utility of the agent. Second, liquidity 
crises are commonly accompanied by major market drops (cf. Russia in 2008, US in 2001, Germany 
in 1945) which lead to significant unavoidable losses for the agent. Third, during a trading inter- 
ruption, agents are unable to liquidate their holdings to finance consumption. Thus, agents may 
experience shortage of funds and are forced to dramatically curtail their consumptions while their 
investments are frozen. This was widely documented in the mass media in autumn of 2008 and 
was one of the major mechanisms of the transmission of the financial crisis into the real economy. 

In this paper we investigate the effect of such trading interruptions on the optimal investment 
and consumption strategies of a power-utility agent. We maintain the Black-Scholes model for the 
risky asset and the assumption of frictionless trading outside of the liquidity event. By keeping our 
setup identical to the classical Merton model except for the feature of liquidity freezes, we are able 
to draw clear comparisons to the original benchmark. 

Because historically major liquidity crises have been rare events, we undertake asymptotic analy- 
sis of the effect of illiquidity on optimal strategies. Indeed, due to their rarity, statistical estimation 
of the frequency of liquidity crises is fraught with difficulties and the corresponding model parameter 
cannot be well calibrated. One of our contributions are closed-form formulas for the first-order cor- 
rections to the classical Merton solutions for power-utility agents. These formulas show the impact 
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of the liquidity crises frequency and duration and provide simple rules of thumb for the investors- 
Simila r asy mptotic expa,nsion s for other models of liquidity were undertaken in [Rogers and Zane 
(12002) and iMatsumotol pOOfil ^ . 

An alternative viewpoint is that liquidity crises are commonplace but have a very short dura- 
tion. Mathematically, this means that liquidity is a fast-scale process, modulating the slow-scale 
geometric Browni an motion i nodel o f the risky asset. Such theory of fast-scale hybrid diffusions 



was developed by 



Il'in et al, 



( 1999 



I). 



The effect of fast-scale stochastic volatility on optimal 



investment-co nsumption was studi ed bv lBauerle and Riederl (|2004l ) in a regime-switching volatility 



model and by iFouque et al.l (j200ol ) in a diffusion setting. The key observation is that the singu- 



lar perturbation of the HJB equations requires averaging of the Sharpe ratios rather than of the 
volatilities themselves. Our results below (fully explicit for the logarithmic and hyperbolic utilities) 
complement this analysis; in particular our perturbation is of a new type since the trading control 
is only active in one regime. On a practical level, this homogenized model allows investors to adjust 
their strategies to account for a fast-scale friction in the market. 

Our model is closest to the recent work of lDiesinger et al.l ()2010l ). In the latter work, which was 



one of the inspirations for our research, the authors consider a model with similar motivations and 
techniques. Following their setup, we treat liquidity as a two-state Markov chain which modulates 



the m arket dynamics and imposes constraints on the trading strategies. While iDiesinger et al 



(j20ld ) work on a finite horizon and study maximizati on of terminal u t ility, w e instead focus on an 
optimal consumption model on infinite horizon. Also, IDiesinger et al.l (j2010l ) do not carry out any 
asymptotic an alysis and do not treat fast-s cale liquidity. An earlier study of illiquidity effects was 



carried out by 



Schwartz and Tebaldi 



(j2006l ) who considered optimal investment where the liquidity 



shock times are known in advance. Finally, our use of a finite-state Markov chain to modulate 
the investment regimes al s o res oi iates with regime-swi tching models of optimal investment, see 
Sotomavor and CadenillasI ((20091) ; IZariphopouloul (|l992l ) . 

The rest of the paper is organized as follows. Section [2] gives a formal description of our hybrid 
diffusion market model and the optimization problem for the investor. Section [3] provides the 
solution for a log- utility maximizing investor. Section H] repeats the analysis for a general investor 
with HARA utility; in particular Section 14.21 provides explicit solutions for the case of hyperbolic 
utility; also implicit formulas are obtained for square-root and inverse-square-root utilities. All these 
analytical results are illustrated with several numerical examples in Section [5l In Section [6] we study 
the fast-scale limit of our illiquidity model and its effect on optimal investment and consumption. 
Section [7] extends to a finite-horizon setting where we consider a log-utility maxi niizing agent, as 



well a s provide the asymptotic analysis of the terminal log- utility maximizer of lDiesinger et al 
(j20ld ). Finally, Section [8] concludes. 
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2. Market with Regime-Switching Liquidity 
Let (Mf) be a two-state continuous Markov chain with infinitesimal generator 

\ ^10 --^10 / 

Associated with M are two Poisson processes A'oi and A'^io with intensity rates of Aqi and Aio 
respectively. The process Ni^i-i, i G {0, 1} can be viewed as a counting process for the number of 
transitions of M between state i and Thus, M has the representation dMt = 1{a/^_=o}^^oi(^)~ 
l{A/,_=i}dA^io(t). 

The process M is a proxy for the market liquidity. When Mt = 0, the market is liquid and the 
usual Black-Scholes model for the risky asset (henceforth called 'stock') price applies: 

(1) dSt = fiStdt + aStdWt, 

with (Wt) a one-dimensional Wiener process on a stochastic basis ($7,7^, P). When = 1, the 
market is illiquid; the stock is non-traded and its price experiences a deterministic exponential 
growth at rate a < r: 

dSt = aSt dt. 

Moreover, at the times of transition from Mt^ = to Mt = 1, the stock price experiences an 
instantaneous drop of L%, St = {I — L)St^. This jump represents the abrupt decrease in stock 
price due to a catastrophic market event (such as 9/11 in the USA, collapse of Lehman Brothers, 
etc.). Overall, the stock obeys a hybrid jump-diffusion model 

(2) dSt = fiMt St dt + GMt St dW{t) - LStl{M,_=o}dNoi (t) , 
with jl = [fi, a]'^ , a = [a, 0]-^. 

2.1. Wealth Dynamics. We consider an economic agent who solves an optimal investment- 
consumption problem in this regime-switching market. Namely, at any given date t the agent 
(i) invests a fraction vrj of her total wealth in the stock; (ii) invests the remainder of her wealth 
in a bank account paying interest at fixed rate r; (iii) consumes at rate ct per unit time. We 
make the additional assumption that consumption is only financed through cash. This introduces 
an additional cash crunch constraint whereby the agent may run out of money to consume during 
an illiquidity shock. 

Denote by (Xj) the total wealth of the agent at date t. Then (Xt) evolves according to 

— = r(l -TTt)dt + TTt — - — dt 
-^t Jt 

(3) =(r--^)dt + l{M, =o}{Mf^-r)dt + 7:tadW{t)-L7:tdNoi{t)} + l{M, =i}7Tt{a - r) dt. 

■^t 

The instantaneous jump of size L in the asset price at the beginning of a liquidity shock has a 
double effect: it decreases overall wealth, and also increases the fraction of money held in cash. 
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mitigating the cash crunch. Overall, if today's wealth is Xt- with a proportion vr^- held in stock, 
an instantaneous downward jump of L% in the stock price causes the total wealth to decrease to 
Xt = ((1 — ■Kt-L)Xt- and the fraction invested in S to shrink to vr^ = g{iTt-) where 

^ 7t{1-L) 

= 

Note that g{-) is an increasing convex function from [0, 1] onto [0, 1] which reduces to the identity 
map q^tt) = TT when L = 0. 

Remark 1. Our model can be straightforwardly extended to the case where the jumps of S at 
transition times of M have an arbitrary discrete distribution /l(-) which is bounded away from 
1 (and such that its realizations are independent of the rest of the model). All the results below 
would continue to hold in that case after replacing with g(vr) = 7r(l — -E'[L]) (1 — 7r£'[L])~^. We 
maintain constant jump sizes for ease of notation and interpretation. 

In contra-distinction from the classical Merton setting, the number of shares irtXt/St held by the 
investor has liquidity constraints. Namely, admissible trading strategies are required to maintain 
the above quantity constant on the time intervals where = 1. This represents market illiquidity 
whereby the agent is completely unable to trade stock, and corresponds to exogenously given 
dynamics of vr^ when Mj = 1. To make sure the agent does not go bankrupt, the constraint 
TTt G [0, 1] is imposed. Thus, no sho rt-sales or leveraged p ositions are allowed. Conditional on 



Mt = 1, the dynamics of (vr^) are (cf. (jPiesinger et al.l . 120101 . Lemma 2.1)) 



geu D C 
201ol . ] 



(4) dTTt = TTt [(1 - vTj ) (a - r) + ct/Xt] dt, 

where the second term represents the increase in proportion of wealth in stock due to consumption 
that is paid for with cash. 

Our aim is to quantify the effect of market illiquidity on optimal consumption/investment. In 
particular, two asymptotic regimes are of interest. First, the case Aqi being small represents a 
marginal possibility of liquidity breakdown. Observe that with Aqi = and Mq = we recover 
the classical Merton problem. Therefore, the asymptotics of the limit Aqi — ?• can be viewed as a 
regular perturbation of this model that takes into account the agent's response to illiquidity. 

Second, a homogenized limiting model is obtained by re-scaling the generator Q of M to be 
Qe A Q_ £qj. g ^ i_ This represents a market that exhibits a fast scale of liquidity regimes that 
impose additional constraints on the trading agent. In other words, both Aqi and Aio are large, 
with a fixed ratio A = representing the proportion of time spent in liquid regime. We also 

rescale the jumps as L(e) = Le which in the limit becomes a negative drift on the stock. In the 
limit e — )• we obtain a singular perturbation of the original Bl ack-Scholes model. Such singularly 



perturbed hybrid diffusions with fast switching were studied by lll'in et al.l (|l999bl ). 



Remark 2. In the classical Merton problem with strictly positive consumption, the only constraint 
on trading is no short-sales, ttj > 0. With liquidity freezes, we have the additional no-leveraging 
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constraint vr^ < 1. Thus, even as Aqi — >■ or e — )• above, if vr > 1 our value functions do not 
converge to the original Merton solutions. Consequently, in markets with large Sharpe ratios, a 
liquidity gap would be present even with negligible chance of market freeze. See Section 13.21 and 
Remark U] below. 



2.2. Consumption on Infinite Horizon. As our main model, we consider an agent maximizing 
cumulative utility of consumption on infinite horizon. We assume that consumption is absolutely 
continuous, so that we can define a consumption rate q with respect to the Lebesgue measure on 
M+. The objective function is then given in terms of time-additive utility of consumption. 



where p > is the inter-temporal substitution factor for consumption and u : IR+ — )• M is the strictly 
concave utility function satisfying the Inada conditions lim^^^o ^'(2;) = +00, u' (x) = 0. In 

([5]), the expectation is with respect to measure conditional on initial endowment Xq = x. Below 
we will focus on the case of Hyperbolic Absolute Risk Aversion (HARA, also known as CRRA or 
constant relative risk aversion) utilities, u{x) = ^ with 7 € (— oo,0) U (0,1) and u{x) = logx 
for 7 = 0. With HARA utilities, investment proportions are known to be independent of current 
wealth which also carries over to our model and simplifies the interpretation of our results. 

Admissible strategies for ([5]) are defined as follows. The observable filtration is always F = 
(J-f) generated jointly by S and M. The investment process (vr^) is required to be an J^-adapted 
process with E[Jq tt^ ds] < 00 for any t < 00. Next, we also require that (ct) be J^-adapted with 
E[J^ e-P^\u{ct)\ dt] < 00. In the case where u{0) = —00, ruin is infinitely costly and therefore any 
admissible strategy must additionally satisfy P(7r( < 1) = F{Xt > 0) = 1 for any t < 00. We denote 
by ^0 the set of all admissible trading strategies that satisfy the above constraints, and finally 
by ^ C ^0 the subset given in ([5]) that in addition satisfies the freeze constraint during liquidity 
shocks expressed in @. 

To solve ([5]), we employ the usual Dynamic Programming paradigm. Denote by V^{x) the value 
function of the agent with initial wealth endowment of x and conditional on Mq = 0; similarly 
denote by V^{'it,x) the value function of an agent with an endowment of x and a fraction vr of 
her wealth in stock, conditional on Mq = 1. Standard stochastic control arguments suggest the 



(5) 



sup 
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following Hamilton-Jacobi-Bellman equation 



sup\-pV^ + {x{r + (a - r)7r) - c)V^ + 7r(l - 7r)(a - r)V^ + vrc/xV;,^ 



c>o 



1 



(6) 



< 



sup 

7rG[0,l],c>0 



+ Aio(y°(x) - VHtt, x)) + u{c)} = 0, 



+ Aoi(Fn 



vr(l-L) 



(l-7rL)x)-y°(x))+'u(c)} 



= 



1 - ttL 



limV (vr, x) = ti(0), and 



limy"(x) = n(0) 



In the illiquid regime 1 only consumption can be optimized; the liquidity shock expires at rate 
Aio and once it is over liquid trading resumes and remaining utility will then be V^{x). In the 
meantime, the fraction invested in stock evolves according to (HD and total wealth Xt according to 
([3]). Analogously, starting in the liquid regime the agent optimizes the positive fraction of wealth 
invested in the stock and consumption. All the while she is anticipating the possibility of liquidity 
shocks that occur at rate Aqi and entail immediate loss of L% in the stock value. The different 
dynamics of (Xt) under the two liquidity regimes lead to modified derivative terms with respect to 
X in ([6]). 

We now state the following verification theorem; the proof is given in the Appendix. 

Proposition!. Suppose thatJ°{x) £ {R+) , {n , x) G C^'^(IR+ x [0, 1)) are two smooth functions 
with polynomial growth in x that satisfy the HJB equation ([6]). Then J^{x) > V^{x), J^(7r, x) > 
y^(7r,x). Moreover, if there exist continuous functions tt* ,c^'* ,c^'* satisfying 



such that the strategy {tt,c) is admissible, then J'^(x) = V^{x) and J^{ir,x) = V^{-k,x) and 
n* ,c^'* ,c^'* are optimal strategies. 

2.3. Asymptotics for Small Aqi- For the asymptotic case Aqi — 0, we may uncouple the two 
value functions above. Indeed, on finite horizon when Aqi is sufficiently small, the probability of 
more than one — )• 1 transition of (Mj) is negligible. In fact, the probability of n such transitions 
is approximately (using Aqi ^ Aio) 0{e~"'^"^) = o(e~'^°i) for n > 1. Thus, assuming that Mq = 0, 
we need to consider at most two liquidity regime-changes in (Mt). For this purpose we will define 
the value functions V and W as follows: V{x) represents the maximum value an agent can extract 
when Mq = and no liquidity shock has occurred so far; W{7r, x) represents the maximum value 
when Mq = 1 so the liquidity shock is already present and after its conclusion we revert to the 
classical Merton setting (with value function V{x)). 




c^'*{x) G argmaxj— cJ^ + u{c)}, c^'*{tt,x) G argmax{— cJ^ + vrc/xJ^ + u{c)} 
c>0 c>0 
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Lemma 1. We have V^{x) = V{x) + XoiV{x) + ©(Aqi) where 

sup <^ -pV + (x(r + {fi- r)7r) - c)K + 7Ta;^vr^o-^T4a; 

7r6[0,l],c>0 ^ 

+ Aoi(VF(^^^i^, (1 - 7rL)x) - + ^(c)| = 0; 

(7) < 1 - vrL J 

sup-^ — pVK + {x{r + (q — r)7r) — c)Wx + 7r(l — 7r)(a — r)W7r + ttc/xWtt 
c>0 *- 

+ Aio(y(x) - Win, x)) + n(c)} = 0. 

Proof. Consider the formal asymptotic expansions V^{x) = V^^{x)+XoiV^^{x)+XI^V^'^{x)+0{Xqi) 
and V^{iT,x) = V^^ (n , x) + XqiV^^ {tt , x) + X^^V^^ {n , x) + O (X^^) . Plugging into dH) and matching 
powers of Aqi shows that V^{x) = V{x) is the solution of the classical Merton problem, while 
V^^{tt,x) = W{Tr,x) solves the second equation in d?]). Matching terms of O(Aoi) in the equation 
for shows that V^^{x) = V{x) solves the first equation in ([7]). □ 



3. Log Utility 



We now specialize t o the case u{x) = logx. We first recall the following classical result (see 
(jKaratzas and Shrevd . 1 19981 . Chapter 4) ) . 



Proposition 2 (Example 9.24, p. 150 in lKaratzas and Shrevd (jl998l )). The solution of the optimal 
consumption problem with log-utility is given by V{x) = ^logx + h with 



(8) 



' r/p — 1 + log p 6^ 
p 2p^ 



(p-r) 



Moreover, the optimal consumption/investment strategies are given by c = px, vr 



fi—r 



In the special case a = r, meaning that the excess return of the stock is zero in the illiquid 
regime, an explicit solution to ([U]) is possible. 

Theorem 1. The solution of ^ for u{x) = log x and a = r is given by (vr, x) = ^ log x + h^ir) 
and V^{x) = ^ log x + b where 



(9) 

and 

(10) 



h{n) 



-/p-l + log p + Aio6 + log (1 - 7r(i+^io/^)) 
P + Aio 



r/p-l + \ogp p + XiQ 

b = h Q[iT ) 



P 



p(/0 + Aoi + Aio)' 



(11) C(^*)^ sup 

7re[0,l] 



{p - r)7r - TT^ay2 X 



+ 



A, 



log (1 - <7(vr)(i+^-/^) ) + ^ log(l -nL)\. 
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The optimal investment strategy in the liquid regime is the maximizer vr* in (jlip and the optimal 
consumption schedules are 

(12) c^'*{tt, x)=px-{l- vr(i+^w/p)). 

(13) c°'*(x) = px. 

Theorem [1] presents a complete solution of the two-state liquidity model for log-utility. The 
two asymptotic regimes will be studied in the next subsection; in the meantime let us make a brief 
remark about (jl2p - (|13p . First, observe that optimal consumption in the liquid regime c^'* is myopic 
in the sense of being unaffected by liquidity shocks. In the illiquid regime, c^'*(-,x) IS 
reduced due to the cash crunch. This curtailment of consumption increases as vrt t 1 fmd the agent 
becomes more liquidity-constrained. The effect decreases as the duration of liquidity shocks shrinks 
(Aio t) or intertemporal discounting decreases {p \.). 

Theorem [T] also explains the behavior of stock investments in our model. When Mt = 0, the 
agent maintains a constant fraction of her wealth in stock, vr* given in (jlip . When a liquidity shock 
occurs at date r, the illiquid fraction jumps to vr,- = (/(vr*) and then starts growing according to 
TTt = Il{t — t) where 11 is the solution of the deterministic differential equation dll{t) = pll{t){l — 
Il{ty^'^^°^P) dt, n(0) = g(vr*) resulting from combining (jH) and (fT2]) . Once the liquidity shock 
expires, wealth is immediately rebalanced to vrt = vr* again. 

Proof. Comparing with the form of the Merton value function in Proposition [2l we make the ansatz 
given in the statement of the Theorem, V^{Tr, x) = ^ log x + /i(7r) and V^{x) = ^ log x + b. 

Plugging-in into ([6]) we find that h must solve the nonlinear first order ode 

(14) -{p + Xio)h + r/p-l + Aio6 - log(p-i - 7r/i'(7r)) = 

with b the constant in and corresponding optimal consumption level of 

c ' (7r,2;) = — — -. 

p — 7r/i'(7r) 

Taking vr = removes the differential term and we find 

r/p-l + log/j-F Aio6 



(15) h{Q) 



P + \ 



10 



However, since zero is a singular point for (jl4p with /i'(0) = 0, the value of /i(0) is not sufficient to 
specify uniquely the solution of (fT4|l . Indeed, the general solution of (fill) is 



(16) hU) = h(0) + log f 1 - Ctt^^+^^'^/p)) , 

P + Aio V / 

with integration constant < C < 1. To find C, we invoke the cash-crunch constraint which implies 
lim^-t-i V^{iT,x) = —DO or lim^^i H^tt) = lim^|i /i'('/r) = — oo. Subject to this restriction we obtain 
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C = 1 in (jl6p which combined with the expression for /i(0) in (jlSp yields ([9]). Furthermore, 
leads to 



h'iTT) 



p(l _ -n-(l+Aio/p)) 

and together with the first-order-condition for c^'* implies (|12p . 

Returning to the liquid regime and substituting the ansatz of V^{x) we obtain that b must solve 

-(p + Aoi)6 + {r/p - 1 + logp) + Aoi/i(0) + C(^*) = 0, 

with (■(•) defined in (jlip . We note that C('^*) is independent of b; on the other hand recall that 
/i(0) in (jlSp is a function of b. Simplifying we end up with (jlOp . 

Since C(^) contains all the terms involving vr, its maximizer vr* must be the optimal investment 
strategy. When Aoi = 0, the maximizer is tt from ([8]). Consequently, because h'{Tr) < for all vr, 
the optimal investment fraction vr* < vr is less than the classical Merton fraction in Proposition [2l 

Finally, the ansatz for implies immediately that c^'* = {V^)~^ = px which matches the 
consumption schedule of Proposition [2j □ 

When a ^ r, the ode (I14p has another h'{Tr) term and is no longer separable. As a result, no 
closed-form solution is possible; however the rest of the conclusions of Theorem [1] continue to hold, 
e.g. (jlOp . (|13p . The condition a = r can be interpreted as zero excess return for the stock in the 
illiquid regime. By adjusting the value of L, the negative jump at the beginning of a liquidity 
shock, the general a ^ r setup can be effectively approximated. 

To interpret Theorem [H we compare the liquid value fu nction V^jx) and the classical Merton 
solution V{x) in terms of the efficiency loss introduced by Rogersl ( 2001 ). Namely, let < < 1 



be the percentage of total wealth that needs to be subtracted from an agent facing the classical 
Merton problem to make her utility equal to that of the agent facing liquidity shocks. Then for 
log-utility, = 0(^°s) satisfies 

-\ogx + b=-log({l-G^^°^^)x)+h ^ 0('°s) = i_exp(p(6-/i)). 



P P 
Simplifying the above expression using (llOp we find 

(17) 0(^-) = 1 - exp + Cin*) / + ) . 

\ 2p p + Aoi + Aio/ 

Thus, the loss in utility depends on the function (^(•) in (|lip . as well as on the squared Sharpe ratio 
6*^ of the stock. 

3.1. Small Aoi Asymptotics. We recall that tt* satisfies the first order condition from (jlip . 

L 



{p-r) - CT^TT* + Aoi 



ph'{g{7r*))9'{7r*) 



1 -7r*L 



0. 
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When vr < 1, we may expand around tt to obtain an asymptotic expression of the optimal 
investment fraction tt* for Aqi small: 



TT* = TT — AoiTTi + 0(Aoi) where 



(18) =Afr^+ 



o-2ll-7rL (l-7rL)2 i _ 5(^)(i+Aio/p) I ' 

In a typical market p ^ Aio and therefore g(7r)'^i''/ is negligible; it follows that the effect of liquidity 
shocks on stock investment mainly depends on L. 

To understand the small Aoi-asymptotics of the value functions we use d?]). This implies that to 
leading order, we should replace b with h in the expression of h^n) in (jlSp and therefore b satisfies 

-pb + {r/p - 1 + log p) + Aoi/i(0) + C(^*) -Xoih + 0(Agi) = 0. 

Plugging- in the expansion of vr* in p8]) and simplifying leads to b = h + Aqi^i + 0(Ao^) where 



P 



r/p- 1 + logp + Aio/i f 1 1 n , ^ wi ^Mi+Aio/pN 

h + - log(l - vrL) H — log(l - g{TT) ^ ^°'^) 



P + Aio p p + Aio 

When p <^ Aio, the above implies that the asymptotic efficiency loss is 

(19) e = Aoi0i + 0(Ag,) with ©i^^-i^Mi^. 

2Aio P 

Thus, the asymptotic efficiency loss ©i increases quadratically in the squared Sharpe ratio of the 
stock and has a logarithmic relationship with respect to stock jumps during liquidity shocks. 

3.2. Assets with Large Sharpe Ratios. The above expansion is valid if tt < 1, i.e. the Sharpe 
ratio of the risky asset is not too large. If tt > 1, then the constraint vr* < 1 remains in force 
throughout and it is the case that lim;s^,)^_j.o ^* = 1- Supposing that g{TT*) = 1 — e for some small 
parameter e and using 1 — (1 — e)" ~ ae for any a > the first order condition for vr* reduces to 

// - r - - — - Aoi— ^ + 0(Aoi) = 0. 

(1 + Aio/p)e(l - L) 1-L 

Therefore, 1 — vr* is asymptotically proportional to Aqi and 



Using this expansion inside V^{x) we find 



(21) V\x) 



1, r/p— 1 + log p a — r — a/2 

- log X H \ 5 

p p p2 



+ 1 , ,\ , Aoi logAoij + 0(Aoi) 
I p(p + Aio) J 



where the term in the square brackets corresponds to the value of a Merton strategy with vr = 1. 
Note that the leading term in the expansion is now of order Aqi log(Aoi) due to the cash crunch. 
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4. HARA Utility of Consumption 

We next consider the case of HARA power utilities, u{x) = where 7 G (—00, 1) \ {0} is the 
risk-aversion coefficient. We start with quahtative analysis of the nonlinear ode satisfied by the 
value function in the illiquid regime and then derive the small Aqi asymptotics for the liquid-regime 
solution. The remaining subsections deal with special values of 7, whereby closed-form solutions 
are possible for W{'tt,x). Throughout this section we make the standing assumption that asset 
value during the liquidity shocks grows at the riskless rate a = r. 



Proposition 3 (Example 9.22, p. 149 in iKaratzas and Shrevd (119981 )). Let u{x) = x^/j, 7 S 
(-00,0) U (0,1). Define 



5 = p — -yr 



2(1-7) 

If J > and 5 < then V{x) = +00 and infinite utility of consumption can be extracted. Otherwise, 
the Merton optimal investment/ consumption problem has the solution 

(22) V{x) = f ■—, where / ^ /^^ ~ 



7 \ ^ 

The optimal strategies are 

23 7r = - cx = - X. 

Below we will work with the power-type utility functions that admit a convenient scaling property 
in wealth. Namely, in analogue to the form of the classical Merton solution, we have 

Lemma 2. For 7 G (— oo,0) U (0, 1), the value functions are homothetic of degree 7 in x. 

Proof This standard result follows from the linearity of the wealth dynamics ([3]) which imply that 
if( cf ) is an admissible consumption schedule for the problem with initial endowment x, then (/3cf ) 
is admissible for initial endowment of f3x, /3 > 0. Moreover, the admissibility set A of trading 
strategies is independent of the endowment. Suppose now that (c^'*) is an e-optimal consumption 
control for V^{x). Then using above. 



V"{f3x) > E 



00 



e-P'iPcf*yds 







> /3TT/°(x) 



e. 



Applying the same inequality with (cf^'*) gives that V^{x) > /? '^V^{jx) — e and since e is arbitrary, 
we must have V^{l3x) = f]^V^{x). Similar arguments work for V^{7r,x). □ 

We first concentrate on the case of small Aqi and work with the uncoupled equations ([7]). Based 
on Lemma [21 we make the guess W{7r,x) = f{TT)x^/j. Simplifying shows that / solves 

(24) i-p + jr)f + (1 - 7)[/ - TTj-'f']^ + Aio[/ - /] = 0. 
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The boundary conditions for (j24p will be discussed shortly. For tt = 0, we can directly find /(O) as 
the solution of /(O) = H~^{0) where 

(25) H{x) ^ i-p + 7r - Aio)x + (1 - j)x^ + Aio/. 

Making the substitution e~^ = tt, (f){z) = (p — 77' + Aio)/(e~^) — Aio/ we transform (fM|) into 

(26) V(^) = (p_^r + Aio)f^V '^^-./.(z)-Aio/, z>0. 



7 Vl-7. 
The last equation is a nonlinear separable first-order ode on with the right-hand-side resembling 
the Bernoulli ode except for the constant term Aio/- Its special form allows for much tractability. 

Consider first the case 7 < 0. It can then be checked that H{x) in (j25p is increasing on [/, 00) 
and H{f ) < so that ()25p has a unique root with /(O) > /. Transferring into ()26p we obtain 

hm (p{z) = {p-^r + Aio)/(0) - Aio/ = (1 - 7)/(0)^, 

and the ode (j26p has a horizontal asymptote. As for log-utility, the cash constraint implies that 
V^{l,x) = —00 which means that /(I) = (/)(0) = -|-oo. Moreover, we have 4>'{Q) = —00 and 4>{z) is 
monotone decreasing on [0,oo). 

Next, consider 7 > 0. When hyperbolic risk aversion is positive we have n(0) = 0, and the agent 
can tolerate zero consumption. Therefore, V^{l,x) and W{l,x) are finite. Taking c = directly 
in the pde ([7]) for W{tt,x) and making the ansatz W{l,x) = /"^(l)^ we obtain 

/+(!) = ^i5^</</(0). 

p--fr + Aio 

Recall that p—^r > is a necessary condition to guarantee finite / and holds thanks to assumptions 
of Proposition [TJ Note that since the marginal utility of consumption at zero is infinite u'{c)\^_^ = 
+00, we still have <j)'{0) = /'(I) = -|-oo. Conversely, the function H{-) in (I25p is now decreasing and 
has a unique root /(O) > /. To summarize, when 7 > 0, </>(•) is finite, with the singular boundary 

condition (p{0) = 0,(j)'{0) = +00, and is monotonically increasing on [0, 00) to the horizontal 

7 

asymptote (1 — 7)/(0)t-i _ 

The following technical lemma, proved in the Appendix, clarifies the structure of ()26p . 

Lemma 3. The ordinary differential equation (I26p has a unique solution satisfying the above bound- 
ary conditions. Moreover, for j < we have that (j){z) (resp. /(vr) ) is of polynomial growth as z ^0 
(resp. vr — 7- Ij. 

Recall that the optimal consumption level is given in terms of (p as 

, ^ <^(-log(7r)) V/^ 
(^,x)= (^^— — J -x. 

By Lemma[3]the consumption level c^'*(7r) is therefore polynomial in (1 — vr). Let ttI'* , with initial 
condition tTq'* = vr^'*, be the illiquid fraction of wealth of stocks when consuming according to 
c^'* and remaining in the illiquid regime for t time units. Then dir]'* = ttI'* c^'* {irl'* , Xt) / Xt dt < 



14 M. LUDKOVSKI AND H. MIN 

(1 — TTt)^^ dt for some Ci. It follows that there exists a constant C2, such that vr^'* < 1 — t~^^ and 
the cash crunch is approached at a polynomial speed. 

To summarize, solutions of ()24p can be written in terms of indefinite integrals. In particular, 
when — 7""*^ is a positive in teger (i.e. 7 = —1, —1/2, —1/3), the ode in (p6|) is of the type known as 



Chini differential equation (jKamkd . 1 19441 ) and may be solvable analytically. Below, we will explore 
further the cases 7 = — 1 (leading to a quadratic ode of hyperbolic type), 7 = —1/2 (leading to a 
variant of the Abel cubic ode) and 7 = 0.5 (leading to a reciprocal ode again of hyperbolic type). 

Returning to the liquid regime, let us denote .^(vr) = / ( t^t—t^ ) " (1 — ttL)"^. We make the ansatz 



V{x) = ^ ■ b which leads to the following algebraic equation for the constant b: 



fo7\ -p + ir - Xoi 1-7 7 J N ^ , 7-1 2 2, , \ ?(^)\ n 

(27) oH OT-i + sup < [fi — rjirb -\ cj vr 6 + Aqi r = 0- 

7 7 ^G[o,i] I 2 7 J 

4.1. Small Aoi Asymptotics. Expanding b and vr* simultaneously in powers of Aqi as 

b = f + Xoibi + 0{Xli) and vr* = vr + Aqivti + 0(Agi), 

we find that for Aqi small the first-order correction to vr* can be written as 

(28) ..=- '"';"";ff +0(Ag.). 

7(1 - f 

Plugging back into (j27p . and simplifying we find 

-S{f + A0161) + (1 - 7)(/ + Aoifti)^ + Aoi(C(vr) - / - A0161) + 0(Agi) = 0, 

so that 

, _ m-f _ m)-m-i) 



(29) 



The derived formula (I29p has several implications. First, translating in terms of efficiency loss 
we have 0^^) = AoiG^^^ + 0(A^i) where 

(30) er^ = — 1 1 - I — ~ 1 ^ '"1 = = f _ (7 - 1) 




7/ V / J 



Again, unless vr is large (e.g. bigger than 70%), ^(vr)// is close to 1 and the correction due to 
illiquidity is small. Also, the consumption during the liquid regime is driven by the same expressions, 

_ ^^^^ , ^ 

/ 

Equation (j30p shows that the liquidity efficiency loss is driven by the key function /(tt) from (|24p . 
In the remainder of this section we investigate special values of 7 whereby closed-form expressions 
for /(vr) (or (t){z)) are possible. 



cO'*(x) = rr • (/ + A0161) '^^^ = c{x) + Aoiy) 
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Remark 3. The coupled equations ([6]) can be treated similarly, by making the ansatz V^{x) = bx'^ /j 
and V^{tt,x) = f{7r)x'^/'y. This leads to the coupled equations for /(vr) and b, analogously to 1^4 
and (l27ll. 



4.2. Hyperbolic Utility. For the case of hyperbolic utility, 7 = —1, the nonlinear ode for the 
illiquid regime ()24p can be solved analytically. 

Theorem 2. Let /? = + r + Aio). The solution of ^ with u{x) = —\ is V^{x) = —B/x, 
V^{t^,x) = —F{'k)/x where 

(31) n-) = ^{1 + ^(-1)^ + • 

(32) (p + r + \m)B - 2Vb + sup |(/i - r)7r5 - a\'^B - -^^^-^(^(vr))) = 0, 

7re[0,l] 1 - ttL J 

and 7?(-l) = Vl + /3AioS. 
T/ie optimal consumption schedules are c"'* = x/B"^ and 

, /3(l-^''(-i)) 
c ' (vr, x) = X- 



^(1 + ^.;(-i))[(l + ^(_l))2 + vr^^(-i)(i _ r^{-l)Y] 

The proof of Theorem [2] follows straightforwardly by checking that the solution -F(vr) satisfies 
the corresponding version of (j24p and then substituting back into the liquid regime equation (j27p . 
The algebraic equation (I32p can be solved easily numerically for B. For small Aqi the solution of 
the uncoupled HJB equation ([7|) is obtained by using r^{-l) = (1 + /3Aio/y/^ in ^ and keeping 
([32]) as is. 



Remark 4. In the case of the large Sharpe ratio tt > 1 we can use (I3ip to re-compute the small 
Aoi-asymptotics. Omitting the details that are very similar to the derivations of ()20p -()2ip we find 

-1/2 



^* = 1 - V^Tvri + o(Aoi), vri = 2^^ {/(/i - r - 2^^)^ 
Note that the correction terms for both vr* and b are now of order \/Aoi- 

4.3. Square- Root Utility. Another closed-form solution for W{tt, x) is possible when 7 = 0.5 
and u{x) = 2-y/x. This case is instructive for showing the interaction of positive risk-aversion with 
the liquidity freezes. 

With 7 = 0.5, the ode (I26p can be re-written in the integral form as 



I dz + C= I „ ^'^ ^- d(l), (/.(O) = 0. 

The integral on the right hand side can be computed using partial fraction decomposition. The 
resulting implicit solution is 

, / 2(f>(z) x-l-(Aio/)/'?(0.5) / 26(z) x-l+(Aio/)M0.5) 

(33) e^= 1 + ^ -A -1 ^ ^ 



r/(0.5) + Aio/^ ^ r/(0.5)-Aio/ 
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where 7/(0.5) = y 2(/9 — 0.5r + Aio) + (Aio/)^- We remark that for typical parameter values, 
limj;^oo 0(-z) is small, so that (j){z) (and hence /(vr)) is very flat and consumption is not very 
sensitive to vr. 

4.4. Connection with Abel Cubic Ode. When 7 = —1/2, the function in the RHS of ([26]) is a 
cubic polynomial with polynomial discriminant D = ^(/3 + 0.5r + Aio) — ^(/5 + 0.5r + Aio)^(Aio/)^. 
In a typical market we have D < and thus this cubic equation has one real root (namely, hi) 
and two complex conjugate roots (namely, /i2 and /i2). These roots can be obtained from the cubic 
formula. Let p = 3f?(/i2) and q = 3'(/i2) > denote the real and imaginary parts of the complex 
root /i2. Then 

7]{-0.5) , 9 __1 _V3/, 9 

'"(p + 0.5r + Aio) 8r7(-0.5)' ^" 2 ^' " 2 4r/(-0.5) 



where 



vi-0.5) + 0.5r + Aio)2/3 ( 4Aio/ + 1 / (4Aio/> ^ 



1/3 



4^^ ' -^"^ \^ -^^^ ^"^^ (p + 0.5r + Aio); 

Since is a separable ode, it can be again solved explicitly using the above roots and partial 
fraction decomposition. Furthermore, since the two complex roots are conjugate, we find the real- 
valued solution to the Abel ode. To conclude, the implicit solution to the ode ([26]) for 7 = —0.5 is 
given by 

(34) lorf 1 + ^^"^ arctanf ^ 

= -^{{hi - pf + q^){p + 0.5r + \w)z. 

Remark 5. For the case of 7 = —1/3, the degree of the polynomial ode (j26p is four. The roots of 
a polynomial up to degree 4 can be explicitly found in terms of its coefficients. Using the quartic 
formula, the polynomial ode (j26p of degree 4 can then also be explicitly solved. 



5. Numerical Examples 

To illustrate our derivations above we present a series of numerical examples. Recall that the 
model for 7 = admits closed- form solutions given in Theorem [1] and for 7 = — 1 in (j3ip - (|32p . 
These formulas are explicit up to finding the maximizer tt* of the liquid regime that can be easily 
obtained through a standard numerical optimization algorithm. Formulas (|33p and (|34p give im- 
plicit solutions for W{-k,x) for 7 = +0.5 and 7 = —0.5 respectively. Inverting these expressions 
numerically and plugging them into (|27p one can again straightforwardly solve the model. Finally, 
for general 7 one can use numerical integration routine on (j56p and combine it with a root-finding 
algorithm for (f27ll to compute V{x) and W{-k^x). The coupled system ([6]) can be solved similarly 
by employing a Picard iteration over successive numerical approximations to V^{x) and y^(7r, x). 
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The key parameter in our model is the coefficient of risk-aversion 7. To ihustrate its effect 
we focus on the three cases of 7 E {0,0.5,-1} which correspond to log-utihty, square-root and 
hyperbohc risk-aversion coefficients. As discussed above, if vr is moderate (less than 60 — 70%), 
the effect of liquidity freezes is typically negligible; therefore the interesting situation is where vr is 
close to 1. Since vr depends on 7 (cf. (|23p ). to compare different HARA utilities we fix vr and then 
adjust the excess return /x — r as 7 varies. In other words we take fj, = r + 7rcj^(l — 7). 

Our base parameter set is p = r = 0.05, a = 0.167, a = 0.05, L = 0, Aoi = 0.1, Aio = 2. The 
stock drift is fi = 0.05 + 0.025(1 — 7), equivalent to tt = 0.9, which is rather large and makes the 
effect of illiquidity noticeable. This base case is supposed to represent a realistic stock market that 
experiences a major liquidity crisis about every ten years (e.g. historically in 1987, 1998, 2008, etc.); 
once started, the crisis lasts for about six months. To study the effect of liquidity shock parameters, 
we also vary the A's and L. 

To study our asymptotic approximations for small Aoi, in Table [T] we compare the solution using 
the coupled system ([6]) to the approximate formulas (fTSj) (7 = case) and (f29l) (7 / case). The 
last two columns also compare the efficiency loss (relative to the classical Merton optimizer with 
-^01 = 0) using the exact solution and the small- Aqi approximation. As expected, when Aqi i 0, 
the optimal investment proportion vr* converges to the Merton solution vr and the quality of the 
asymptotic approximations improves. Same conclusions hold for the efficiency loss. 

As can be seen in Table [H for realistic parameter values and L = 0, the efficiency loss is mild 
(less than 2% for logarithmic and hyperbolic utilities and less than 1% for square-root utility). 
This conclusion is consistent with the " relaxed investor" property of the Merton problem found by 



Rogera (|200ll ) ; iRogers and Zand (|2002l ) who studied efficiency loss when trading is discretized on 



a time grid. Thus, inability to rebalance continuously is not crucial, while cash crunch constraints 
will also be weak unless Aio is small. 

Not surprisingly, if liquidity shocks are accompanied by significant price decreases, a major utility 
loss is experienced since the excess return of the asset is effectively cut. Even if a liquidity shock 
occurs once every 10 years and causes a relatively minor negative jump of L = 10% in asset price, 
the optimal investment fraction vr* in the risky asset is reduced from 90% to 52% for log-utility and 
69% for hyperbolic utility (the square-root investor is barely affected and down- weighs to 88.8%). 

In general, the derived asymptotic expressions underestimate the adjustment in risky investments 
TT* < vr -|- AoiVTi and overestimate the efficiency loss < Aqi©!- This is intuitive given the 
convexity /concavity of the respective quantities and the linear approximations taken. Note that 
comparison of given quantities for different utilities is rather difficult, although we can claim that 
the square-root optimizer is highly insensitive to liquidity and the hyperbolic optimizer is somewhat 
more sensitive than the log- utility agent. Overall, Table [T] seems to indicate that while the risk of 
negative asset jumps is crucial, the risk of illiquidity per se is not a significant driver of optimal 
investment decisions. 
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Log Utility 

Case 


u{x) = log{x), fl = 
vr* vr + Aqivti 


0.075 

Aoi0i 


Base 

w/ Aoi = 0.05 
w/ Aoi = 0.02 

W/ Alo — 4 
w/ L = 0.1 

w/ Aoi = 0.5,Aio = 10 


0.879 0.846 

0.886 0.873 

0.893 0.889 
n nnn n qoo 

u.yuu u.oyy 
0.521 0.467 
0.900 0.900 


1.08 1.15 
0.56 0.58 
0.23 0.23 
U.04 U.OO 
13.97 18.13 
1.06 1.11 


Hyperbolic U1 

Case 


ility u{x) = —x^^, 
vr* vr + AoiVTi 


= 0.1 

Aoi0i 


Base 

w/ Aoi = 0.05 
w/ Aoi = 0.02 

W/ Alo — 4 
w/ L = 0.1 

w/ Aoi = 0.5,Aio = 10 


0.857 0.683 
0.867 0.792 
0.879 0.857 
u.oyo u.yui 
0.691 0.652 
0.900 0.916 


1.92 2.36 
1.01 1.18 

0.43 0.47 

n on n qq 
u.yu u.yz 

17.16 18.55 

1.77 1.83 


Square-Root Ut 

Case 


ility u{x) = 2y^, /i 
vr* 71" + AoiVTi 


= 0.0625 

Aoi0i 


Base 

w/ Aoi = 0.05 
w/ Aoi = 0.02 
w/ Aio = 4 
w/ L = 0.1 

w/ Aoi = 0.5,Aio = 10 


0.895 0.898 
0.897 0.899 
0.899 0.900 
0.898 0.900 
0.888 0.335 
0.878 0.838 


0.589 0.623 
0.302 0.312 
0.120 0.125 
0.304 0.314 
16.737 20.213 
0.606 0.638 



Table 1. Optimal investment proportions and efficiency loss for HARA utility- 
maximizers with 7 = (log- utility) , 7 = —1 (hyperbolic utility) and 7 = 0.5 (square- 
root utility). measures the percent efficiency loss compared to Aoi = (Merton 
problem) and 0i is the asymptotic efficiency loss derived in p!7|) and (pOj) . 

To further explore the role of investor's risk aversion, Figure [U compares the optimal consumption 
schedules in the illiquid regime for three different levels of 7, namely 7 G {0.5,0,-0.5}. While 
comparison of value functions is hard due to different utility units, all consumption is proportional 
to total current wealth x and moreover must decrease as the proportion of wealth in cash (1 — vr) 
shrinks. Recall that in the classical Merton setting, c* /x > p when 7 < and c* /x > p when 7 > 
and the proportion of wealth consumed is maximized at 7 = —1. 

As shown in Figured! c^'* is decreasing in 7 for moderate values of vr and is increasing in 7 for 
TT close to 1. This confirms our intuition that the cash crunch is most felt by the most risk-averse 
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7r 



Figure 1. Consumption schedules c '* in the ihiquid regime. Parameter values are 
fixed at a = 0.166,^ = 0.0625, r = p = 0.05 and Aio = 1,L = 0. We compare 
the square-root (7 = 0.5), logarithmic (7 = 0) and inverse square-root (7 = —0.5) 
maximizers. All values are normalized by current wealth (i.e. we plot c^'*{tt,x)/x). 

investors. Interestingly, around vr ~ 0.9, no obvious ordering between c^'*(x;7) is possible. For 
vr < 0.8, /'(vr) is very small all 7's, and therefore consumption is essentially constant in vr. 

6. Homogenized Limit 

We now consider the homogenized limit e — t- with respect to the transition matrix = Q/e 
of M. This corresponds to the model where liquidity regimes change on a fast time-scale, with 
^ = '^io/(-^oi + -^10) proportion of time spent in the liquid unconstrained regime. Formally, we 
assume that that the transition rates of M are \io/e and Aoi/e, the jumps are rescaled as L = eZ, 
and then proceed to carry out a first-order formal asymptotic expansion in the small parameter e. 

As e — 7- 0, the effective duration of liquidity shocks goes to zero. Therefore, the investment propor- 
tion process vrt converges to a constant value, which is the optimizer of the equation defining V^{x). 
Simultaneously, this means that V^{-,x) becomes less and less sensitive to vr, lim^^QV^{'iT,x) = 
for all TT G [0, 1) (super-polynomially), and therefore the optimization over vr in ([6]) is simplified. 

As e — )• 0, the ergodic Markov chain M averages out the drift and volatility of S. Let us 
call a cycle a pair of transitions of M from state to state 1 and back to state zero. Then the 
length of a cycle is about ~ ^{-^oi + ^ki^ ) during which the expected mean of an S'-increment is 
~ H^Xqi +aeXiQ — eL and the quadratic variation of S over one cycle is ~ eXQ^a^. Simplifying the 
above expressions which become exact as e — )• and using the fact that the homogenized process 
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remains a diffusion we obtain that the averaged asset dynamics S are 

dSt/St = {a + {fi-a- XoiL)X) dt + X^'^a dWt. 

However, as we wih see below, the optimal strategy will not converge to the Merton solution 
corresponding to S. 

6.1. Log Utility. First, consider the log- utility investor. We directly work with ^ and make the 
ansatzy°(x) = -Bo+eBi+0(e2) + l log x and ^^(vr, x) = BQ+eF{TT)+0{e^). The leading term is the 
same in both equations because the dominant term in the expression for is Xi^iy^ {x) — V^ {i: , x)). 
Collecting terms of order 0(1) in the illiquid regime we obtain 

-pBo + -(r + (a - r)7r - p) + \ogp = Xio{F{7t) - Bi). 
P 

Back in the liquid regime, we find vr^^^ = where 

Ghom = -(/^ -r- XoiL + ^(a - r)) 
a Aio 

and Bq solves 

(35) i?o = ^/^-^ + ^°g^ +A%. 

Comparing with the benchmark in ([8]), we may interpret the limiting value function in (I35p as 
solution of the classical problem for a risky asset with a modified squared Sharpe ratio of XOf^^^ 
rather than 9^ as in ([8]). Therefore, the homogenized loss in utility computed in (|17p translates into 



(36) lim ^ 0jf = 1 - exp(^(AC^ - 9^)). 



Counter-intuitively, we find that the derivative of 0/,'°^] with respect to Aqi and Aio may be of either 
sign depending on parameter values. For instance, \i p> a> r and L = then d@^^^/dXiQ > 0; 
however if p is large and a <r then dQ^j^^f^/dXiQ < 0. 

In the special case a = r, the explicit results of Proposition [T] apply and yield ([5^ directly. 
In that case since the growth rate of the risky asset during liquidity shocks is a = r, we can 
assign the illiquid regime a Sharpe ratio of zero. Then X9f^^^^ corresponds exactly to averaging 
the squared Sharpe-ratios of each liquidity regime according to the invariant distribution of M. In 
contrast, the averaged price process S has excess return of A(^ — r — Aoi-Z>) and volatility of \fXa. 
Thus, the illiquidity constraint causes the agent to apply the wrong strategy over the averaged 
stock process. This is particularly clear in the case L = whereby lim^^Q vr* = vr which is the 
optimal strategy of the liquid regime that is employed over the homogenized price process S. 
In the related setting of a fast-scale reg ime-switching stock model with M-modulated drift and 



volatility (but no liquidity constraints), iBauerle and Riedeij (j2004l ) showed that the log- investor 



also loses 0^'°^^ of utility compared to the Merton model. Therefore, for the log-investor in the 
homogenized limit the liquidity constraint completely disappears. When a 7^ r, 9hom no longer has 
an averaging interpretation and the liquidity constraint remains in force even in the homogenized 
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limit. For Example 1 in Section [5] above, we find Q^^'^f^ = 1.066%, which is very close to the value 
e = 1.061% for Aoi = 0.5, Aio = 10 in Tabled) 



6.2. HARA Utility. To obtain the homogenized limit for the HARA investor, we carry out a 
similar formal asymptotic expansion of equations (|24p and (j27p in the small parameter e. Starting 
once more from the ansatz V^{x) = {Bq + eBi + 0{e'^))x'^ /j and V^{7r, x) = {Bq + eF{7r))x'^ we 
obtain 



t 7 — 1 



Aio(F(7r) - Bi) = i-p + 7r + (q - r)Tr)Bo + (1 - 7)^0 > 
and the optimizer vr* is 



1 



vr 



Denoting 6 



(7) A 1 
hom 



(1 - 7)^2 



^ if! — r — XoiL + ^^(a — f)) , the equation for Bo reduces to 



fi-r - XoiL + 



A, 



01 



Aio7 



(a — r) 



5n = 0, 



Bo = {1 - j)'-' { p - jr - XieZ'J 



(7) n2 7 



7-1 



2(1-7), 

Again, we interpret the expression for Bo as the solution of the classical Merton problem after 

9(7) 

hom' 



modifying the Sharpe ratio of the risky asset to A^^^^. The resulting efficiency loss is 



(37) 



^hom ^ 



I p-jr- X^^^^ 



For Example I in Section[5]the efficiency losses for the homogenized limit are = 1.742% for 



Hyperbolic utility (compare with 1.766% loss for Aqi = 0.1, Aio = 2 in Tabled]) and = 0.600% 

for Square-root utility (compare with 0.606% loss in last row of Tabled! for 7 = 0.5). Thus, these 
limiting expressions give accurate approximations even for moderate values of Aqi and Aiq. 



7. Logarithmic Consumption on Finite Horizon 

We proceed to study a finite horizon version of the model above. With finite horizon, time 
inhomogeneity introduces additional effects. In particular, close to terminal date T, both the cash 
crunch and the opportunity cost of illiquidity vanish. The corresponding stochastic control problem 
under consideration is now: 



(38) 



sup 



e log Cs ds 
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Proposition 4. For model (j38p . the Merton solution is 
(39) V{t,x) = k{t) log x + h{t) 

= ^(1 - e-^(^-*)) log X - ^(1 - e-^(^-*)) log 



(/i — r)^ r — p 



1 - 



T/ie optimal consumption and investment strategies are: 

px 

^ " 1 - e-p(^-*) ' 



TT 



fi — r 



l _ e-p(^-*)' 



i + p(r-t)) 



Assuming that Aqi is small and with a slight abuse of notation, denote by W{t, ir, x) and V{t, x) 

the first-order correction terms to V^{t,TT,x) and V^{t,'7r,x) as in ([7]). Then the finite-horizon 

analogue of ([7]) is now 

f c 
sup-^ Wt + {x{r + {a - r)7r) — c)Wx + tt{1 - 7r)(a - r)WT, + tt-W-,, 

OO X 



+ Xw{V{t, x) - W{t, vr, x)) + log c - pW^ = 0, 
sup <j Vt + {x{r + {p- r)7r) - c)T4 + ^x^vr^cj^Fxi 



7rg[0,l],c>0 



+ Xoi(.W{t,^^^—^,{l-7rL)x)-V{t,x)) + logc-pV]=0, 



with zero terminal conditions W{T, Tr,x) = V{T,x) = for all vr, and the singular boundary 
condition W{t, l,x) = — oo for all t,x. 

Theorem 3. The value functions for (I38p are V{t,x) = k{t)logx + h^{t) and W{t,'7r,x) = 
k{t) log X + h^(t,7r) where /i* satisfy 

(40) = hl- ph^ + rk{t) - 1 + XioCHt) - h\t,'K)) - log{k{t) - 7r/i^(t, vr)), 



(41) h"{t) 



~ip+Xoi){s-t) 



[rk{s) - 1 - log k{s) + ^(s)] ds, 



with terminal conditions h^{T) = 0, h^{iT,T) = and 



(42) m = sup {fi - r)7rk{t) - -a'n'kit) + Aqi • 

7rg[0,l] ^ 

The small Xoi-asymptotics are 



k{t)log{l-nL) + h'{t,^^^—^) 



(43) 



TT = TT -|- 



^01 



^ +±-hl{t,g{n))g'{n) 



1 - ttL fc(t) 

T 



h^{t) = h{t) + A, 



01 



log(l - vrL) 



(/ii(s,(7(vr))-/i(s))ds 
'l_e-/'(^-*)(l + p(r-t))' 
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Let us note that the key equation (j40p is a first-order nonhnear pde on [0,T] x [0, 1]. Because 
of the cash-crunch, it has the degenerate boundary condition lim7r->i vr) = — oo which makes 
its numerical solution nontrivial. To overcome this difficulty, in Section 17.11 below we carry out an 
asymptotic expansion of h^{t, •) around vr = 1. 

The small Aoi-asymptotics of the optimal investment fraction show that the effect of liquidity 
shocks mainly depends on stock jumps since is negligible for typical market parameter values 
and large time horizon. The correction term in the small Aoi-asymptotics of the value function h^, 
however, depends on stock jumps as well as efficiency loss during liquidity shocks. This correction 
term is always negative and decreases in magnitude as t — )■ T. 



7.1. Cash Crunch on Finite Horizon. When the market is illiquid and vr ~ 1, the agent 
experiences a cash crunch. The resulting singular boundary condition V^{t, l,x) = — oo presents a 
numerical challenge for solving (j40p . To resolve this issue, we carry out an asymptotic analysis of 
the optimal behavior when (1 — vr) is small. 

Lemma 4. Let h^{t,7r) be given in (|30|) . Then asymptotically for small 1 — vr and any < d < 1 

1 _ Q-{p+^io){T-t) 

(44) h\t,l-d7r) = h\t,l-7T) + ^ log d + o{l - tt) . 

P + Aio 

Proof. When most of the wealth is in the stock and Mf = 1, consumption is driven by the amount 
of cash available, y = {1 — Tr)x. Indeed, because the agent must maintain positive consumption, 
she concentrates on avoiding a complete cash depletion. Let r ~ Exp{Xiq) be the next jump-time 
of M. Then the agent solves the optimization sub-problem 



W{t, 1 - vr, x) = C{y; x) = supE 

ict) 



fTh{T~t) 

/ e-''^ log Cs ds + e-''("^(^-*)V°(r, + Y^) 

Jo 



where (It) is the amount of cash at date t obeying dYt = (rYj — q) dt, and Xt is the wealth in stock 
obeying Xt = xe^*. Since (1 — vr) is small, Xr ^> Yj- and therefore consumption is asymptotically 
linear in remaining cash, since Cg effectively only appears in the logc^-term. Therefore, for any 
fraction d < 1, 

C{d • y; x) ~ / e"''" log(d • c^e-^^''' + Aioe-(''+^")"y°(s, xe") ds 
Jo 

= C{y;x)+ [ {e~^P^^'°> log d}ds 
Jo 

= C{y; x) H — log d. 

P + Mo 

Re-arranging the last expression in terms of we obtain (j44p . □ 
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Figure 2. Relative utility of W{t, vr, x) with respect to x) for different propor- 
tions of cash holdings 1 — vr. On the vertical axis we plot the percent efficiency loss 
= 1 - exp . Here T = 2,r = p = 0.05, ^ = 0.075, = 1/6 and L = 0. 



We may use ()44p as a numerical boundary condition to solve (I40p . Namely, we construct a finite 
grid {0, Att, 2A7r, . . . , 1 — Avr} and solve (140 p using a finite-difference scheme, applying 

1 _ Q~{p+^io)(T-t) 

h^{t, 1 - Avr) = h^{t, 1 - 2 Att) log 2, 

P + Mo 

which is now a finite boundary condition at vr = 1 — Avr. Such a numerical solution is presented in 
Figure [2] which compares the resulting vr) to the Merton solution h{t). As t — )• 0, the percent 
efficiency loss is linear in vr (becoming proportional in available cash); as t — )• oo, the efficiency loss 
converges to the infinite- horizon formula (jl7p . 



7.2. T erminal Utility vi^itliout Consumption. A related model recently appeared in lDiesinger et al 
(|201ol ). 3KS hereafter. The latter paper was one of the inspirations for our analysis and considered 



optimization of terminal wealth under liquidity shocks. In contrast to our analysis, DKS did not 
allow for consumption but did allow a more general form of multiple liquidity shocks (rather than 
a two-state Markov chain). Similar to our results, DKS obtained explicit solutions for the case of 
log-utility. In this section we carry out the asymptotic analysis of their model which provides an 
instructive counterpart to our optimal consumption problem. 

The DKS objective is 
(45) sup E*'^ [log(XT(l - l|M^=i}vrT^))] • 

Thus, if the market is illiquid at terminal date T, total wealth is reduced by L%. 
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Proposition 5. The Merton solution for maximizing terminal wealth is 
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(46) 



sup E*'-[log(XT)] = logx + (r + - 



Denoteby J°(t,x), Ji(t, TT, x) the value functions corresponding to (|45p . Then ,P solve the system 
sup |x(r + (^-r)7r)J° + ^xW2jO^ + Aoi(J^(t,x,7r)- J°(t,x))l =0, 

7r6[0,l] I ^ J 

j/ + x(r + (q - r)7r)4 + Aio(J°(t, x) - J^{t, vr, x)) + 7r(l - 7r)(a - r) J^''' = 0, 
j'^(T, x) = log X, and J^(r, vr, x) = log x(l — vrL). 



When a = r, the solution reads (jPiesinger et al.l . l20ld . Section 5) J^{t,x) = logx + /"(t) and 
(t, vr, x) = log X + /-"^ (t, vr) where 



(47) 



r + (/i - r)vr2)^^c; - ;Jo-2(vr2)^^c;)2 + ^(?^ + Aio log(l - tt*dksL)) 
I Aio 



1 _ e-(^oi+Aio)(T-t) 



Aoi + Aio 



Jl^g- Aio (T-t) 

Alo 



-(Aoi+Aio)(T-t)> 



A, 



01 



Aoi + Aio 



+ 



Apir 

La?, 



(-1 + e 



-Alo(T-^)^ 



10 



\io 



Aoi + Aio 



!\t, vr) = j\xiofis) + (r + Aio log(l - vrL)))e-^^«(^-*) + log(l - vrL)e-^io 
with optimal investment proportion ir'^xg defined through 



(T-i) 



(48) 



2_* 



= (/U — r) — (7 vr 



DKS ~ Aqi 



1 - Lvr 



DKS 



We proceed to understand the asymptotics of the above expressions when Aqi is small. First, 
expanding ([I8|) in Aqi we find (subject to the condition 1 > {fi — r)/a'^ > 0) 



T^DKS 



a'^ + {^l- r)L - ^((j2 - (/i - r)L)2 + 4^2 L2 A 



01 



2ct2L 



(49) 



fi — r 



A, 



L 



01- 



O(Agi) 



cr2 — (/i — r)L 

As expected, risk of illiquidity reduces holdings of stock, with the effect increasing as (i) a decreases 
or (ii) increases or (iii) the liquidity jump loss L increases. For typical parameter values the 
denominator above is quite small and therefore the relative changes in vrj^^^^ can be quite big even 
for a small change in the underlying parameter. Note that if L = then the optimal strategy 
is always equal to the Merton solution vr, due to the form of the first-order condition of t^*£)xs- 
This phenomenon of myopic investment is unique to log- utility. In the same vein, vr|)^^ is in fact 
independent of Aio to first order. 
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Comparing the expression (j49p of the optimal investment proportion to a similar expression (jlSp 
derived in Section [3] we see the impact of consumption on investment. Indeed, the two problems are 
identical except that intermediate consumption introduces the cash-crunch feature which increases 
the impact of illiquidity. Analytically, this difference is exactly represented by 

T^DKS - ^ = Aoi - 



As already discussed, unless vr is very close to 1 or Aio is small, this difference is relatively small 
and therefore there is little difference between maximization of terminal wealth and maximization 
of utility of consumption. 

Plugging into the expression for in ()47p and simplifying we end up with 

\2 



fit) = {r + ^^^)(T - t) + Aoi • (l - e"^-(^"*) - X^T - t)) 



+ log(l-^^^^)(r-t) +0(A2 



a2 



01; 



Comparing with the limiting case Aqi = in ()46p . we obtain that the percentage loss in utility is 
0(DKS)(^) = Aoi0S''''''^(t) + 0(A2J with 

(50) Sr'\t) = • (e-^-(^-) + A,o(r - t) - l) - log (l - (T - t) 

Both terms above are always positive and decrease (approximately linearly) in magnitude as t — t- T. 
Again, the squared Sharpe ratio makes an appearance. Counterintuitively, as Aio — oo, the 
illiquidity correction does not vanish, even if the amount of time spent in the illiquid regime 
is negligible. This is because there is always a possibility of an illiquid shock right before T, 
proportional to Aqi which would produce a loss of L% of equity. 

Similar analysis can be done for the fast-scale liquidity freezes: 
Lemma 5. Rescale the generator of AI as Q'^ = e~^Q. Then as e ^ we have 



vr* ^ (^-^-/"^-"^ fit) ^ (T - t)r + A(r - ,) (^i^l^^ 
Therefore the homogenized percent efficiency loss is 

(51) &iZ'\t) = 1 - exp ({T - t)^{-xe^ - e^} 



2 



We omit the proof of Lemma [5] which is similar to the computations for (j36p and follows from 
expressions (I48p - (147p after tedious arithmetic simplifications. Also note that the efficiency loss in 
([STI) essentially matches (p6]l except for the discount factor. 
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8. Conclusion 

In this paper we considered a perturbed version of the classical Merton optimal investment- 
consumption model. Our modification accounts for possibility of liquidity shocks during which 
the agent cannot trade and cannot use her stock holdings for consumption. The liquidity shocks 
also lead to negative jumps in risky asset price and reduce its excess return. We obtained explicit 
formulas for HARA utility-maximizers, especially for the special cases of 7 = 0, ±0.5,— 1. The 
economic interpretation of our results was carried out by deriving the percent efficiency loss relative 
to the respective Merton problem solution. We showed that the asymptotic effect of illiquidity can 
be captured by appropriately modifying the Sharpe ratio of the underlying asset inside the classical 
formulas. This gives a simple rule to the investor in terms of accounting for liquidity shocks. 
Overall, we found that for realistic parameter values, the efficiency loss is primarily driven by the 
jump parameter L. 

Our model can be easily extended if one is willing to sacrifice analytical tractability. Natural 
extensions would allow for more complex stock dynamics (such as general jump-diffusion processes, 
non-zero volatility in the illiquid regime), arbitrary utility functions, further market frictions or 
multiple stocks. The resulting HJB equat ions would then be modifie d version of the respective 



classical optimal investment models, see e.g. lKaratzas and Shrevd (|l998l ). Also, other risk objectives 



might be considered to judge the attitude of the investor towards illiquidity. Our first setting where 
probability of a liquidity shock is small resembles the lifetime investment problem, in which case 
quantities such as probability of ruin become important risk metrics. Another direction that we 
leave to future work is to consider the effect of liquidity shocks on indifference pricing of contingent 
claims. 
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Proof of Proposition [T] 

Proof. The proof is l argely standard and follo ws the classical Verification Theorem for stochastic 
control of diffusions (jFleming and Soneri . Il993l . Theorem IV. 5.1). To simplify the presentation, we 
consider the case u{0) > — oo. In the sequel we will primarily work with HARA utilities u{x) = x'^/j 
which are infinite at the origin for 7 < 0. That case can be treated similarly to the proof below by 
using the resulting horn othetic property of the value fun ction (see Lemma [2] and the met hods in 
Shreve and Soneii (j 19941 )). Let us also mention the paper ISotomavor and CadenillasI (j2009l ) which 
had a similar verification theorem for optimal consumption in a regime-switching model. 

Define u{x) = u{x) — n(0). Then u is another utility function, and K[f^ e~P^u{ct) dt] = 
E[Jq°° e~''*n(ct) (it] — u{0)/p. Moreover the admissible strategy sets are the same for the prob- 
lem with u(-) and u{-). Therefore, without loss of generality we assume for the remainder of the 
proof that n(0) = 0. 



To apply the conclusions of ([Fleming and Soneri . 1 19931 . Theorem IV.5.1), the following steps need 
to be verified. 

Step 1: Check that V^{-) < 00 are finite. Case 1: sup^u{c) = +00. In this situation of "positive 
risk aversion" , it may well be the case that the value functions are infinite since even a modest level 
of consumption can generate a lot of utility; see Proposition [3] below. To avoid this we assume that 
the liquidity-unconstrained problem with value functions 



sup E 



n(ct(Xf'^''))(it|Mo = 



and similarly for V^{7r,x), are finite. We refer to ISotomavor and CadenillasI ( 20091 ) for exhaustive 
analysis of that setting. Since A Q Ao we immediately obtain that under the above assumption, 
V^{x) < V^{x) < 00 and V^{tt,x) < V^{Tr,x) < 00. 

Case 2: sup^u{c) < 00. In that case, we immediately have V^{7r,x) < V^{x) < {supcn(c)}p-^ 

Step 2: We next show that J*(-) solving ([6|) are upper bounds for the value functions. 

Let (vr, c) be an admissible strategy. The assumptions on (vr, c) then guarant ee that the sde in ([3l) 
has a unique strong solution (in fact, a closed-form representation is possible, see lSotomavor and Cadenillas 

(booa)). 

The basic idea is now to apply Ito's lemma to the semi-martingale {X^'^''^, Mt). Let 6 = inf{t : 
-^x,n,c _ Since wealth is required to stay non-negative, no consumption is possible on {6,00] 
and it follows that V^{0) = V^{7r,0) = 0. Fix initial wealth x and constants 0<ai<x<a2<oo 
and 03 < 1. Let r = ti A T2 A T3 A 6*, with 



n ^ inf {t : X^'^'' < a,} , T2 ^ inf {t : '"'^ > aa} , ^ inf {t : {nt > 03} n {Mt = 1}} . 
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For any t < oo, the quantities X^'"^'"^, J^{Xs'^''^), J^{Xs'^''^), J^{7Ts, Xg'^'^), being solutions of 
([6]), are local martingales. Moreover, they are bounded on [0,t A r] and consequently are all true 
martingales until tAr. Applying the optional sampling theorem we obtain 



(52) 



J°(x) -E 



tAr 



> E 



Now the right hand side in (|5'2p is non-negative, so that J'^(x) > E^[Jq ^ e~''*u(cs) ds] and taking 
fli —7-0,02 — ^ 00)^3 — ^ 1 and using Monotone Convergence Theorem (MCT), we obtain J^{x) > 
W[JI^^^ e~''^u{cs) ds]. Finally, taking t — )• oo and again using MCT we find J^{x) > y'^(x; (vr, c)) 
as desired. Similar steps apply to J^{7r,x). 

Step 3: Verify existence of optimal controls. In our case sufficient first order conditions for the 
consumption optimizer above imply that 

Assumed smoothness of guarantees existence of c*'*. The investment optimizer solves 

ifi - r)xV^ + 7T*x'a^V^, + {g{7T*), (1 - 7t*L)x) g'{7T*) - LxV^ (^(vr*), (1 - 7t*L)x) = 0, 

and an interior solution on [0, 1] is guaranteed if ^ > r. Indeed, in the latter case the above 
expression is positive when vr = and goes to — oo as vr — ?• 1 (we have u'{0) = +oo and as proportion 
of liquid wealth shrinks, consumption must be reduced to zero, implying that lim7r-|-i = — oo). 

□ 



Proof of Lemma [3] 



Proof . The ode ()26p is separable and therefore the corresponding theory may be applied (j Walter 



19981 . p. 31). We begin by re- writing in terms of indefinite integrals as 



(53) dz + C 



7(p -7r + Aio) [jTzryj ) - l(t> - 7A10/ 



First consider the case 7 > 0. An easy argument shows that (/>(•) is increasing and therefore the 
solution satisfying (p{0) = is given implicitly by 

(54) z=\ H{x) = {{p-jr + X^o)i- - x - Aio/ . 



/o -/H{x) [ VI -7 

7 ~ 7 

Recall the horizontal asymptote limz_!,oo (/'(•z) = (1 — 7)/(0) . Checking that ff((l — 7)7(0) t-^ ) = 

Hif{0)) = 0, and using H{x) > for x < /(O) and 7 > we see that lim.z^^i^\ = +00 

and thus the boundary conditions are matched correctly. 
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When 7 < 0, (p{z) is decreasing; therefore for N > (p{oo) the unique solution to (j26p satisfying 
(j){0) = N is given by 

(55) = 

Moreover, since the continuous function H is bounded on the compact domain [(/)(oo), A^], the 
solution (!){•) is Lipschitz on [0, oo). Since the leading power in H{x) is 1 — > 1 for negative 7, 
the integral in (I55p is proper and converges as iV — )• 00. Therefore, taking the limit, we find that 
the unique solution with lim^-j-o 0(z) = +00 is 

-dx 



(56) 



z) jH{x) 

1-1/7 



Fix N' large enough. Then for x > N', jH(x) = j{p — jr + Aio) ( + o{x^~^''^) and 
therefore 

r ^ r W^'-'dx = -^{Nr^\ A^ip-jr + Aio)(l - t)^""' + o{N'). 

Jn' ^H{x) Jn' ^ ^ 

Comparing with ()56p we find that 

0(e) = e^(l - 7)'-^(/0 -ir + Aio)^ + 0(e). 

Thus, (j) grows polynomially (at rate 7) around zero. Finally, transforming back into the 7r-variable, 
and using that for e small enough, e"*^ > 1 — 2e we obtain 

/(l_,) = 0(e7) 

showing that / also grows polynomially as the cash crunch is approached. □ 

Appendix A. Proof of Theorem [3] 

Proof. Once again following the form of the unperturbed problem, we posit that under liquidity 

shocks, the value functions are V{t,x) = {t) log x + h^(t), W{t,TT,x) = (t) log x + {1,11). 

Substituting the ansatz for the illiquid regime value function W we obtain that c^'* = 

and after some simplifications that k^ = k. Moreover, the wealth-independent term h^{t,7r) solves 

(BQD. 

Returning to the liquid regime, the ansatz implies that c^'* = x/k^{t) and 

k^t logx + /i° + rk^(t) - 1 + Aoi logx[A:^(t) - k^{t)] - Aoi/i°(t) 

+ logx-logA;°(t) -/3[A:°(t) log x + /iO(t)] + C(t, ^*) = 0. 

where 

e(t,vr*) = sup |(/i - r)T:k\t) - \a\^k\t) + Aqi [k\t)log{l - ttL) + {t, gin))]} . 

7rG[0,ll ^ ^ ^ 



32 M. LUDKOVSKI AND H. MIN 

The terminal condition is k^{T) = 0, hP{T) = 0. Considering the terms multiplying logx, we find 
again = k. Thus optimal consumption in the liquid regime cP'* is again unaffected by liquidity 
shocks. Similarly, collecting terms of 0(1), we see that solves the first-order ode 

h^t-ph^{t) - \oih^{t)+rk^{t) - l-\ogk^{t)+C{t,TT*) = 0. 

Since this first-order ode is linear and ^(i,7r*) is independent of h^, we can apply variation of 
constants to obtain (j4ip and (j42p . Next, from equation (j42p we know the optimal investment 
fraction tt* satisfies 

(57) if^-r)- a\* + Aqi 

Using equation ()57p we expand vr* to the leading order for small Aqi to obtain the first half of (I43p . 
Plugging-in this asymptotic expression for vr* into and expanding (|4ip for small Aqi to the leading 
order, we recover the second half of (j43p . □ 
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